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THE  EFFECT  OF  THE  SURFACE  LAYER  OF 


ABSTRACT: 


A  STAR  UPON  BLACK  BODY  RADIATION 
EMERGING  FROM  THE  INTERIOR 


Murray  Turoff 
Brandeis  University 

A  form  of  Bunnann's  Theorem  is  used  to  obtain  an  expansion 
of  the  radiation  intensity  at  the  surface  of  a  star  as  an 
e3q)ansion  explicitly  dependent  on  the  black  body  function, 
absorption  coefficient,  density,  and  the  first  few  deriva¬ 
tives  of  these  functions  with  respect  to  depth*  Conditions 
for  the  convergence  of  the  series  are  established*  The 
e^qpansion  is  modified  for  the  case  of  organized  outward 
motion  of  the  surface  layer  and  the  effect  of  this  motion 
on  an  absorption  line  is  examined*  Conditions  for  the  shift 
of  the  line  "center'*  to  lower  or  higher  frequencies  due  to 
the  effects  of  the  orgauiized  motion  and  the  decay  of  teipera- 
ture  in  the  surface  layer  are  also  examined* 


Introduction 


In  the  outer  layers  of  a  star  the  marked  decrease  in  the  density  and 
temperature  leads  to  the  formation  of  absorption  lines*  The  actual  shape  of 
the  absorption  lines  we  observe  is  complicated  considerably  by  such  effects 
in  the  outer  layers  as  presstire  broadening,  Doppler  broadening,  the  presence 
of  a  temperature  gradient,  arul  organized  motion*  Although  these  effects  can 
be  taken  into  account  through  a  suitable  choice  of  the  absorption  coefficient 
and  soTirce  function,  the  resulting  integral  esqpression  for  the  radiation  in¬ 
tensity  is  often  difficult  to  treat  analytically.  In  PART  I  of  this  paper 
we  obtain  a  more  convenient  form  of  the  radiation  intensity  by  use  of  a 
lemma  employing  the  property  of  marked  decrease  in  the  density  udiich  occurs 
in  the  outer  layers*  Itfe  also  establish  general  conditions  for  the  conver¬ 
gence  of  the  resulting  eiq)anslon*  In  PART  II  we  treat  a  specific  model 
which  allows  us  to  re-establish  our  conditions  of  convergence  in  terms  of 
the  commonly  used  scale-heights  of  density,  temperature,  and  pressure* 

In  PART  ni  we  modify  our  original  expression  to  take  into  account 
the  presence  of  organized  outward  motion*  The  shift  and  shape  of  the  line 
is  examined  in  detail  and  an  esqpression  for  the  shift  of  the  "center”  of 
the  line  is  derived*  In  PART  IV  the  method  employed  to  eoqpand  the  radia¬ 
tion  intensity  is  used  to  obtain  a  similar  expansion  for  the  radiation  flxix* 
The  lemma  which  is  used  several  times  in  this  tpproach  is  derived 


in  the  Appendix 
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PART  I;  Expansion  of  the  Radiative  Intensity 


Consider  a  stellar  atmosphere  under  the  following  assumptions: 

(l)  The  surface  or  outer  layer  of  the  star  begins  at  x  ■  0  and  extends 
to  X  -  oo  •  The  density  in  this  surface  layer  exhibits  exponential 
decay  as  its  dominant  beha:rior»  i«e*: 


fW  = 

where  essentially  the  scale  height  of  the  density  and 

T  »  -r  if*  for  all  X. 

(2)  The  interior  of  une  star  is  an  approximate  black-body;  therefore  the 
radiation  emerging  from  the  interior  and  incident  on  the  outer  layer 
at  X  -  0  is  given  by  the  Planck  (black-body)  functlon:^'^**^ 


I 

-  I 


(2) 


where  the  temperature  T^  -  T(x"0). 

(3)  The  outer  layer  of  the  star  is  considered  to  be  in  approximate  local 
thermodynamic  equilibrium;  however,  the  decrease  in  both  temperature 
and  density  in  this  region  of  the  star  leads  to  the  presence  of 
absorption  lines  in  the  spectrum  of  the  star* 

Under  these  assumptions  we  may  write  the  radiation  intensity  emerging  from 
the  star  (plane-parallel  atmosphere  approximation) 


1*H=  a 

where 


(3) 


(1.) 

X 

is  the  optical  depth  of  the  surface  layer  and  K(x)  is  the  mass-absorption 
coefficient*  The  first  term  of  (3)  is  the  contribution  to  the  luiiinoslty 
from  the  interior*  ■  cos  9  where  9  is  the  inclination  to  the  outward 


h 


nornal*  Both  K  and  B  are  functions  of  the  frequency  (P)« 

Because  of  the  assumed  exponential  decay  of  the  density,  ve  nay 
enploy  the  relation  (All)  derived  In  the  Appendix  to  expand  (li)  and  obtalnt 

where  for  any  function  of  x,  l*e.  A(x),  ve  define 


Av  = 


^  ±JA 


A  7x 


A  =  -!-  ''V\ 
*'  A  7F 


and 


(6) 


(7) 


The  function  A^(x)  has  units  of  length  and  may  be  thotight  of  as  a  ' quasi 
scale  height* 

We  may  also  apply  (All)  to  (3)  to  obtalnt 

I-  !<£>v 
A 


lf>(i 


where 


(8) 


= 


- 1 


(9) 


+  f(y  -f 

and  the  right  hand  side  of  (8)  Is  evaluated  at  x  ■  0.  We  nobet 

>.^=  (10) 
'  - 

If  T(x)  Is  a  decreasing  function  of  x  as  Is  esqpeeted  In  the  photosphere  and 
lower  chromosphere,  we  note  that: 


Aj  <  A,  »S  S,<  o 


(u) 
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Combining  the  expansions  of  1*  and  ,  we  obtain,  in  lowest  order  t 

A' 

The  second  order  terms  will  be  small  if  condition  (11)  holds  and  the  follow¬ 
ing  conditions  are  met: 


(Hi.; 

«  1 

(13) 

t  (  K,. 

-  •<;)<<  1 

(111) 

T 

(15) 

i  (  Bxx 

(16) 
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PART  II;  Analysis  of  Convergence  for  a  Specific  Model 


To  examine  the  convergence  of  this  expansion  In  more  detail,  ve 
assume  the  absorption  coefficient  Is  the  sum  of  a  continuous  part  and  a 
line  parts^*^ 


with 


K(x)  =  K*{*)  t 

K-M  - 


(17) 


(18) 


L  1^  Q 

where  p  -  - j  y  and  k  are  constants  and  the  factor  (1  -C"*^)  is  the 

HT  c  c 

jP 

correction  for  stimulated  emission*  For  K  (x)  we  adopt  the  following  forms 
which  depend  on  the  part  of  the  line  shape  we  are  considering 


for 


for 


C 


>  h. 


(19) 


(20) 


where  V  ^  is  the  central  frequency,  ^  is  the  natural  line  width,  c  Is  the 
velocity  of  light,  and  k  .  is  a  constant*  $«(x)  is  the  pressure  broaden- 
ing  width,  and  V^(x)  -  thermal  velocity*  Equations 

(19)  and  (20)  give  K^(x)  for  the  center  and  the  wings  of  the  line  respec¬ 
tively  and  are  approximations  to  the  more  general  convolution  integral  which 
is  used  to  combine  both  the  Doppler  and  Pressure  Broadenii^  effects  *^*^ 

We  further  assume 

T(x)=  T.  c' 
f(*)  = 


(21) 
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where  associated  scale  heights  for  these  quantities* 

While  the  assumption  of  a  pure  ejqionential  behanrior  for  the  above  quantities 
may  only  approximate  a  physically  realizable  situation,  it  does,  however, 
allow  us  to  obtain  simplified  conditions  for  establishing  the  validity  of  our 
approach  in  terms  of  the  scale  heights*  We  list  a  number  of  relations  which 
will  be  of  use: 


If 

~  if,  ~ 

(22) 

8,  = 

-  P 

XT(i-erP) 

(23) 

X 

J-  ficP 

(2b) 

K*'  = 

(25) 

K'‘- 

(26) 

i*Sf 

Incorporating  all  these  results,  we  obtain,  after  some  algebra. 


A.  = 


where 


I  *  Xf,  H(r>) 


(27) 


for 


and 


H(v)  = 


for 


_! _ /  jlV;  ?|*  _ 


c 
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In  order  to  satisfy  (13) ,  we  wish  X  ^  H(  p  )  to  be  small  when  H(  D  )  is 

r 

negative*  A  sufficient  condition  for  this  is  that  is  small 

I  pop. 

increases  from  zero  this  condition  is  relaxed.  Vfe  therefore  require  the 
following  inequalities  to  be  satisfied: 


( *<  f  ^  ^  <  1 


(28) 


(29) 


(30) 


(31) 


(32) 
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f(^)‘  sy 


Vfe  see  that  condition  (28)  Is  sufficient  to  make  the  quantities  (32)  to  (3Ii) 
snail.  Vis  also  note  that  (15)  Is  Identically  zero  for  a  pure  exponential 
decay  of  the  density.  We  now  require  the  following  quantity  Is  to  be  snail 


conpared  to  one: 


(Uiff 


It  remains  only  to  point  out  that: 


The  first  two  terns  are  small  because  of  (29)  and  (36);  using  this  fact,  we 


nay  write  the  last  tern  as: 


{K*~  k;  I"*- 


Fron  (2U),  (25),  and  (26),  we  see  this  quantity  Is  also  small  under  condi¬ 
tions  (29)  and  (36).  Furthermore,  It  Is  also  small  when  1. 

either  very  small  or  very  large  conpared  to  unity. 

We  nay  now  sunur.arlze  that  sufficient,  but  not  necessary,  conditions 
for  the  expansion  to  converge  for  siLl  are: 


<<  1 


>>/ 
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The  first  condition  tends  to  restrict  the  result  to  the  forward  direction 
of  the  radiation  (|i  -  l)j  otherwise  it  is  only  a  statement  that  the  outer 
layer  of  the  star  is  optically  thin  to  the  black-bo{^  radiation  ftrom  the 
interior*  The  second  statement  requires  the  density  to  decrease  faster  than 
the  temperature  in  the  outer  layer,  which  is  to  be  expected  for  most  stellar 
photospheres*  The  third  condition  certainly  holds  when  either  one  of  the 
quantities  ^  (natural  line  width)  or  (pressure  broadening  line  width) 
greatly  exceeds  the  other*  Because  of  the  second  condition,  we  can  expect 
that  1  X<  greatest  value  expression  (3)  can  therefore 

obtain  is  approximately  l/8* 

Therefore,  to  a  good  approximation,  under  conditions  (38)  we  may 
use  equation  (12)  for  I^(^)*  Expanding  equation  (12)  for  p.  -  1,  we  havet 

!+(«)?  e(T.)(,-(Cp(A, -A,) j  09) 

We  may  <ilao  define  a  deviation,  ,  from  the  black-body  distribu¬ 


tion  t 


“iw—  '  -'''(‘'■M 


(fit  } 

using  relations  (7)  and  (9)  for  2* 


s  (liO) 
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PART  nil  Effect  of  Organized  Motion 


Ve  now  wish  to  modify  the  original  expansion  for  (8)  to  take  Into 

account  organized  outward  motion  of  matter  In  the  surface  la^rsr*  The  only 

change  necessary  Is  to  replace  B  (  x,V  )  and  K(  x,  P  )  in  (8)  with 

B  (x,  >?(x))  and  where 

P  P 


and 


=  V  ( I  ♦  ^  “(x;] 


(U) 


xO 


u(x)  is  the  velocity  of  the  organized  outward  motion  and  we  define  the  boun¬ 
dary  of  the  surface  layer  (x  ■  0)  as  the  point  irtiere  u  ■  0  and  ^  0» 

For  the  steady  state^  the  conservation  of  current  Implies  pu  - 
constant;  this  can  be  made  consistent  with  the  assumptions  on  the  density 
variation  over  any  finite  interval  by  choosing  u  large  enough*  In  fact, 
since  all  physical  effects  take  place  In  the  first  few  scale  heights, 


violations  of  ^u  >  constant  are  not  important  at  distances  greater  than 
a  few  times  • 

We  may  e3q)and  Sp  snd  obtain: 

^  if!  +  .  ,  , 

From  the  boundary  conditions  we  have  established  In  (1:1),  we  note: 


(1>2) 


(li3) 
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Aasvutdng  ■—<<  1,  we  note  that  in  equation  (8)  for  I*  wo  need  only  modify 
the  first  order  terms;  a  modification  of  the  second  order  terms  for  the 
Doppler  effect  would  give  third  order  corrections.  Vfe  further  see  that  only 
the  derivatives  of  6  and  K  are  modified  which  in  turn  implies  that  only  A  ^ 
and  A  2  modified.  Defining! 


we  have 


^  i-il  tf 

8  •'i/ 

I  - 


JL  JH 
k  77 


e  -  LJv. 

<-77 


A 


Ue  may  then  rewrite  our  first  order  eapression  for  I  (12)  as 


m 


(1*5) 

(1*6) 


b{t^)  e" 


(1*7) 


I  ^  _ 

Af/- 

The  Doppler  modification  of  the  first  order  terms  for  I  ^equation 
(12) J  is  of  the  same  order  of  magnitude  as  the  second  order  terms  for 
found  in  equation  (0),  and  omitted  in  our  ejqpresslon  (li7)  for  I^»  However, 
for  the  purpose  of  attempting  to  isolate  the  effect  of  the  outward  motion 
from  the  other  second  order  terms,  we  will  neglect  this  error  at  present. 

A  quantity  which  is  correct  to  the  second  order  is  the  deviation  (^jj) 
of  the  Doppler  shifted  intensity  (I^)  from  the  normal  Intensity  (I*) .  By 
taking  equation  (8)  both  in  the  normal  form  and  modified  to  second  order 
for  the  Doppler  effect  and  expanding  the  exponentials  and  denominators,  we 
obtain  to  second  order: 
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r  -  -1» 

c 

If  we  assume  that  B  and  K  are  constant  with  respect  to  frequency 


m 


mroi*  f.Via 


width  of  the  spectral  llnej  we  flndt 


at  the  center 


r,  ^ 

on  the  wings.  We  also  find  for  ^  itself,  from  equation  (40): 


(50) 


•  *  J- 


where 


r,. 


(51) 


In  the  diagram  we  have  Indicated  roughly  the  behavior  of  and  1^2 

over  the  line  width.  There  is  no  scale  indicated  as  this  may  vary  greatly 
with  the  respect  choice  of  many  different  parameters.  Although  we  have 
assumed  B  constant  over  the  line  we  see  from  (23)  that  the  factor  - 


By 

decreases  with  Increasing  frequency.  This  has  the  tendency  to  shift  all  the 
deviations  various  amounts  to  the  blue  as  V  increases.  This  shift  is 
partly  due  to  the  temperature  gradient,  since  it  arises  partly  from  the  term 
in  we  see  that  it  acts  opposite  to  the  organized  Doppler  shift.  The 


diagram  Is  only  meant  to  show  that  if  the  Doppler  effect  were  dominant>  l*e. 

-i-ii'  »  _L 
^  rfx  At 

the  maximum  of  the  line  would  be  shifted  to  a  lower  frequency  ( V  <  )  •  If 

0 

the  decrease  in  temperature  is  dominant*  i*e* 

tr  ^  ^ 

we  would  e3?)ect  a  shift  to  the  blue  ( •  It  is  possible,  under  the  right 
conditions,  to  obtain  a  double  peak  in  the  absorption  line* 

In  order  to  try  and  obtain  an  equation  for  the  shift  of  the  absorp¬ 
tion  line,  we  may  expand  (1|7)  under  the  conditions  (38)  to  obtain: 

b(i- -Xl(K,*Bt)) 

f or  -  1  (52) 

¥b  may  also  e^qpand  ^  ^  and  ^  ^  under  the  same  conditions  to  obtain 


>-f(i  *XeKx) 


±1 


(,  .  ) 
-AfBxV  >-\B,  I 


Using  (23)  and  (18),  we  note  that  for^«l[wB  have  approximately: 


Carrying  only  the  first  order  term  yields: 


i; 


B(i-rK) 


($$) 


Tti^rnttire 

Shift 


Doppltr 

Shift 


Potslblt 

Coi^sitt 
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where 


(56) 


Taking  the  derivative  of  (B5)  with  reapect  to  ,  we  obtain  the  critical 
points  in  the  line  at  the  zeroes  of  this  derivative: 


I'  Kr 

For  the  center  of  the  line  we  have,  from  equation  (19) 

77-  ' 

where 


when  ^  <  <  £_  vc  find 


\4e  see  from  (56)  and  (60)  that  under  the  conditipnt 


(57) 


(58) 


(59) 


(60) 


(61) 
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a  shift  to  increasing  frequency  vlll  occur*  If  the  condition  Is  violatedi 
we  will  have  a  shift  to  lower  ^equency*  Ife  note  that  since 

are  both  small  numbers  under  conditions  (38)  we  would  expect  that  — - 

A-  J/i  ^ 

j  •  Therefore  we  may  summarize 


Is  much  larger  than 


or 


that  the  following  conditions  hold: 


j-ia 
«  7x 


J-Jv 

^7^ 


for  a  shift  to  the  blue 


for  a  shift  to  the  red 


(62) 


for  p  small. 

These  conditions  reflect  our  earlier  considerations*  made  on  a  more 
intuitive  basis.  The  singularity  when  -  ^  In  (60)  Is  spurious,  for 

c  </;r  At 

the  expansion  breaks  down  In  that  region. 

For  p  large  we  find  that  A  2  which  means  that  the  second 

term  on  the  right  hand  side  of  (52)  will  be  an  order  of  magnitude  larger  than 
the  third  term.  Therefore  we  expect  the  Doppler  effect  to  be  much  smaller 
for  large  p. 


We  note  that  K 


derivative  proportional  to 


(20)  Is  proportional  to 
-3 


(w 


and  Its 


;  therefore,  we  expect  that  a  similar 


analysis,  using  the  pressure  broadening  coefficient  of  absorption,  leads  to 

j  7+ 

a  third  degree  equation  for  ^  -  0  in 


the  first  approximation.  Us 
therefore  eiqpect  that  the  effect  of  pressure  broadening  Increases  the  likeli¬ 


hood  of  having  a  double  peak  in  the  absorption  line.  This  Is  due  to  the 
extension  of  the  wings  of  the  line  by  the  pressure  broadening.  However, 
equation  (55)  Is  not  accurate  enough  to  examine  this  possibility. 
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PART  lYt  Expression  for  the  Radiative  Flux 


To  Illustrate  further  the  usefulness  of  the  approach  used  here,  we 

derive  a  single  expression  for  the  flux  at  the  surface  of  the  star.  Reeal> 
,1,2,3,U 

I 

a.  .  .1 

(63) 


ling:  ■ 


where  (*0  )  Is  the  net  outward  flux,  we  find  that 


eO 


where 


so 


Ehiploylng  the  lemma  derived  In  the  Appendix  plus  the  expansion  for 
also  stated  In  the  Appendix,  we  find: 


1  +  zKf>  j^Kp^i  +  . .  , 

vhere  the  next  terns  In  the  series  are  of  order  ikpAi  )  $  p 

(Kp>i )  ( Aa.  )•  Equation  (65)  Is  a  relatively  simple  expression  for  the 
flux,  which  may  be  examined  In  much  the  same  manner  as  we  have  treated  the 
Intensity  of  radiation. 
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SuiHBary 

The  method  discussed  here  treats  the  effect  of  the  surface  layer  of 
the  star  on  the  black -body  radiation  emerging  from  the  interior*  If  the 
density  in  the  surface  layer  falls  of  rapidly^  the  expansion  obtained  for 
the  intensity  of  the  radiation  converges*  The  advantage  of  the  expansion  is 
that  it  involves  no  integrals  and  does  not  seriously  limit  the  form  of  the 
absorption  coefficient*  The  equation  for  the  intensity  (8)  may  therefore 
be  applied  to  the  particular  model  or  situation  one  wishes  to  treat*  Fur- 
thermore^  the  results  of  the  expansion  are  esqpressed  as  a  function  of  'x' 
rather  than  the  optical  depth*  The  equation  for  the  intensity  or  the  flux 
as  a  function  of  'x*  is  perhaps  more  applicable  to  problems  governing  tiie 
general  dynamics  of  the  surface  layer  of  the  star  where  the  radiative  prob¬ 
lem  is  only  one  part. 

The  results  we  have  obtained  for  the  shift  of  the  absorption  line 
in  the  case  of  organized  motion  is,  of  course,  limited  both  in  accuracy 
and  the  range  of  applicability*  However,  accurate  results  for  particular 
lines  in  particular  model  atmospheres  can  be  obtained  quite  easily  by  using 
our  basic  equation  (8),  modified  by  (li5)  and  (U6)  for  the  Doppler  correction, 
and  a  good  computer* 
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APPEMDIZt 


k  Special  Fori  of  BurwMin'a  Theorew 


life  wish  to  derive  a  form  of  Bumann's  theorem  which  Is  useful  In 
problems  of  the  type  considered  In  this  paper* 

Given  an  Integral  function  of  x:  ^ 

f  9(xpx' 

X. 

where  g(x')  Is  a  monotone  function  from  the  point  x^  to  the  point  x  >  x^, 
the  following  relation  may  be  derived: 


«  i. 

\  I«i(^)(<?(x)-9(x.)) 

y, 

where  Q(x)  Is  an  arbitrary  monotone  function  In  the  Interval  x^  to  x  and  the 
a^  (x^)  may  be  e3q)ressed  In  terms  of  Q(x),  g(x),  and  their  derivatives^ 
evaluated  at  x^*  Defining: 


(A2) 


and  taking  the  derivative  of  both  sides  of  (Al)  with  respect  to  x,  va  have: 

^(x)  f(x)  =  L  ±2—  ( q(x)  -  (q(x.))  («) 

Letting  X  x^  In  (A3),  we  obtain: 

*  a/*-)  <"•> 

Nov  define 
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Note  the  identity  as  x  We  take  the  derivative  of  both  sides  of  (A5) 

with  respect  to  x;  using  (A2)  we  obtain 


or  equivalently 


(A6) 


The  relations  (A6)  or  (A7)  provide  a  straightforward  procedure  for  successive¬ 
ly  determining  the  a^*s«  The  only  other  stipulation  necessary  is  that  a^<  «o 
as  X  ->  x^  which  in  t\irn  depends  on  the  choice  of  Q(x) . 

If  g(x)  is  exactly  integrable  a  suitable  choice  of  Q(x)  is  equivalent 
to  doing  the  integration,  as  the  series  will  terminate  after  the  required 
terms  are  obtained.  A  detailed  discussion  of  the  convergence  properties  of 
this  type  of  series  may  be  found  in  reference  (6) . 

As  an  example,  consider  the  exponential  integrals: 

Let 


we  obtain 


y*h-l 

for  n  •  1  using 


(x.n-i)*’ 


^  ~  J 


C(ir*n-i)1 

Q(t)  =  }(t) 
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we  find 


^  fi.-i - 1— 


The  first  three  terms  of  these  series  expresses  E^(x)  to  an  accuracy  of  two 
places  or  better  over  the  idiole  range  of  x  from  zero  to  Infinity. 

For  otir  purposes  we  are  Interested  In  Integrals  of  the  type: 


(A8) 


where 


Lint 


(A9) 


Defining 


-- 

9  -- -iX 


%  5  LS: 


A  -- 


we  have 


9  M 


(AlO) 


(All) 


In  this  case  we  chose  Q(x)  -  g(x).  Ve  see  that  If  g(x)  uniquely  equalled  a 
decreasing  exponential  the  expansion  would  terminate  after  the  first  term. 
Physically  we  may  think  of  ^  as  a  scale  height.  Even  though  g(x)  may  be  a 
very  complex  function.  If  the  dominant  behavior  was  that  of  an  exponential  we 
expect  the  series  to  converge  rather  rapidly. 
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